Mesoscopic Spin-Hall Effect in 2D electron systems with smooth boundaries 
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Spin-Hall efltect in ballistic 2D electron gas with Rashba-type spin-orbit coupling and smooth edge 
confinement is studied. We predict that the interplay of semiclassical electron motion and quantum 
dynamics of spins leads to several distinct features in spin density along the edge that originate 
from accumulation of turning points from many classical trajectories. Strong peak is found near 
a point of the vanishing of electron Fermi velocity in the lower spin-split subband. It is followed 
by a strip of negative spin density that extends until the crossing of the local Fermi energy with 
the degeneracy point where the two spin subbands intersect. Beyond this crossing there is a wide 
region of a smooth positive spin density. The total amount of spin accumulated in each of these 
features exceeds greatly the net spin across the entire edge. The features become more pronounced 
for shallower boundary potentials, controlled by gating in typical experimental setups. 

PACS numbers: 73.23.-b, 72.25.-b 



Introduction. Spin-Hall effect manifested in the 
boundary spin polarization when electric current flows 
through a system with significant spin-orbit interaction, 
was recently observed in both 3D 0,11, 0] and 2D systems 
[H| ■ Two mechanisms that lead to the eff'ect are typically 
distinguished. The extrinsic mechanism is dominant in 
3D semiconductors and originates from scattering off im- 
purities d, 0, H, Q • Intrinsic mechanism [13, U | of the 
band-structure induced spin precession can be realized in 
ballistic (disorder-free) 2D systems. 

The intrinsic spin-Hall mechanism is of particular ap- 
peal. However, in two-dimensional electron systems with 
spin-orbit coupling linear in momentum (typical for n- 
doped heterostructures) any scattering that leads to a 
stationary electric current via deceleration of electrons 
by impurities, phonons, etc., will negate the precession 
due to external electric field and result in the exact can- 
cellatio n [l2 . isl 14, 15, HI of the bulk spin-current in a 
dc case [l7|. 

One possible way to avoid this cancellation is to move 
into the ac domain with frequencies exceeding the inverse 
spin relaxation time [l^]. Another possibility is to use 
dc fields but make a system sufficiently small and clean 
(ballistic) so that the electron mean free time exceeds 
the time of flight across the systems. The correspond- 
ing scenario became known as the mesoscopic spin-Hall 
effect 22]. While initial theories of spin-Hall effect in 
infinite systems had addressed such auxiliary quantity as 
spin current (for a review see Refs. [ll,[3]), finite geom- 
etry calls for calculations of spin polarization, a directly 
measurable quantity [20| . 

It is important to emphasize that the edge spin po- 
larization in ballistic systems appears not as a result of 
electric field-driven acceleration of electrons and associ- 
ated with it precession of spins. Instead, it originates 
from their precession in the course of electron motion in 
the boundary potential that provides lateral confinement. 
When populations of left- and right-moving states are dif- 



ferent (due to the applied bias) the net effect of this pre- 
cession results in a non-zero spin polarization. Such edge 
polarization was considered, mostl y by numerical meth- 
ods, in several earlier publications 2l|, [22, l23|, |2J, 
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FIG. 1; A) Profile of a boundary potential U[x). B) Ge- 
ometry of the system: two-dimensional electron gas {x > 0) 
is attached to two ideal reflectionless metallic leads filled by 
equilibrium electrons up to difi'erent chemical potentials. Two 
types of boundary scattering are shown schematically. C) 
Spin-orbit-split subbands structure. D) Difference in popu- 
lation of left- and right-moving electrons due to the applied 
bias eV . 

The total amount of accumulated spin near an edge 
of a 2DEG is independent of the boundary potential, as 
found in our recent study [2^ . yet only second order in 
the spin-orbit coupling constant. In the present paper we 
show how to overcome this limitation and achieve a much 
larger local spin-Hall polarization over extended strips 
along the edge of 2DEG provided that smooth bound- 
aries are utilized. Increasing the width of the edge by 
making the boundary potential progressively smoother, 
one can increase the amount of spin accumulated within 
each strip. The sign of the polarization in these strips 
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will alternate so that the net spin accumulation will be 
in agreement with Ref . jl^ . 

Spin dynamics in smooth potentials. Consider ballis- 
tic two-dimensional electron gas (2DEG) attached to two 
ideal perfectly conducting leads which are kept under a 
voltage bias V, see Fig.[T] Hamiltonian of the system fea- 
tures spin-orbit coupling of the usual "Rashba" type 27| , 



value of the force exerted by the boundary potential, 
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We assume that the spin-orbit coupling strength A is 
the same inside 2DEG and in the leads (or, equivalently, 
that switching-on of A happens adiabatically as electrons 
travel from the leads towards 2DEG). In this case the 
applied bias transforms into the difference of Fermi ener- 
gies for the left and right moving electrons far from the 
edges of 2DEG, Ep^ - Ep^ = eV. 

The density of the out-of-plane spin Sz can be related 
to the dynamics of the in-plane spin via simple identities. 



■^^^Syij = -'llj^[H,Sy]lp- 



Vjy , [H,Sy]^2iXszPy. (2) 



Here the operator of spin current = -^[{VipY Syip — 

'0^Vsj,'0] + x'^^'^^a^' ^^"^ ~ i~^y^ ^x)- Since the mo- 
mentum along the edge py is conserved, and we are in- 
terested in the stationary spin-density, we find 



J — ( 



^\{x')Sz^^{x')dx' = ^ — — 



(3) 

Here the sum is extended over all occupied states. The 
exact value of the out-of-plane spin is thus determined by 
the local current J^ix)- Next, we notice that in a smooth 
potential U{x) the subband index is almost conserved. 
The electron's spin remains largely in-plane, while being 
also perpendicular to the momentum, with the local value 
of Px{x) determined from the energy conservation. One 
can now introduce the nonequilibrium (Fig. ID) electron 
distribution for the local value of the boundary potential 
U{x) and the Fermi velocity vpix) — yJ—2U{x)/m and 
calculate the expectation value of J|(a;). This yields. 
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Knowledge of the local distribution, however, does not 
allow one to describe quantum oscillations due to the in- 
terference of waves with different px contributing to the 
same eigenfunction tpi. Thus Eq. @ describes only a 
smooth semiclassical part of spin accumulation, as de- 
noted by the bar, S. We note that the accumulated spin 
S{x) is a convenient quantity for numerical calculations, 
since the integration over x reduces the relative weight of 
the quantum oscillations. Differentiating Eq. ^ over x 
we find that the spin density is proportional to the local 
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{2Tr)^m{vp{x) — X^)vp{x) dx 



(5) 



Note that the net spin polarization across the edge is 
independent of the shape of the boundary potential, 
S{oo) = —X'^eV/12Tr^Vp, and is expressed via the bulk 
value of the Fermi velocity vf{oo). Here by a; = oo we 
assume a point deep inside the 2DEG but yet far from 
its opposite edge. The latter has spin accumulation of 
the same absolute value and opposite sign. 
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FIG. 2: Integrated spin density S{x) — s^dx for the 
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potential U{x) 

curves for a = 8, 4, 2,1 x 10"'^ are shown in yellow, green, 
blue and red respectively. The horizontal coordinate is scaled 
differently for different curves, as x is measured in units of 
lO^a X h/mX. Narrow black lines stand for the semiclassical 
prediction, Eq. Q. The logarithimc behavior, ~ logo, of the 
dip at at U{x) — —mX^/2 {x — 500) is clearly seen. Inset 
magnifies the region near the edge of 2DEG {x ~ 0). 

Numerical solution. To illustrate the predicted depen- 
dence (|4][5]) we performed exact numerical calculation of 
the spin density accumulated near the boundary approx- 
imated by a linear potential U{x) = —Fx with the con- 
stant force F. The smoothness of the boundary implies 
that F ^ m^X^ /h. Fig. [2] demonstrates an excellent 
agreement between the predicted integrated spin density 
Eq. ^ and exact numerical simulations for different val- 
ues of the slope of the boundary potential. 

According to Eq. ^ we find two regions of different 
smooth spin behavior. First, within the strip where < 
vf{x) < X spin density is negative (which is seen in a 
downward slope of the integrated density S{x) in Fig. [2]). 
Farther away, Sz{x) changes sign for vf{x) > A, where 
both Sz (x) and S{x) decrease gradually with increasing x. 

The most interesting is the behavior of spin at the bor- 
ders of these regions, vf = and vf = X. At vf{x) = X 
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the accumulated spin S{x) Eq. ([4]) diverges logarithmi- 
cally. This singularity originates from the accumulation 
of classical turning points taking place when the conical 
crossing point in the spectrum of the Hamiltonian p]). 
see Fig. IC, passes through the Fermi energy. This sin- 
gularity is regularized as 5 ^ logF, according to Fig. [51 

Yet more peculiar is the behavior of both Sz{x) and 
S{x) at the edge of 2DEG, near the point where vf{x) = 
0. The smooth part of the accumulated spin, Eq. ([4]), 
has an infinite jump here (from Eq. ^ it follows that 
5(-|-0) = oo, while of course S{—0) = 0). Development 
of such jump with decreasing slope of the potential is 
seen in the inset in Fig. [21 The jump in S{x) corresponds 
to the formation of a narrow strip with extremely large 
values of spin > along the border. This behavior 
will now be analyzed in more detail. 

Semiclassical analysis. Classical dynamics of electrons 
with Rashba spin-orbit interaction is described by the 
Hamilton function [29l 
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for the two spin-split subbands. The boundary poten- 
tial in Eq. ([6]) is again approximated by the linear func- 
tion. The family of classical trajectories generated by the 
Hamilton function Eq. ([6]) , shown in Fig. [3l demonstrate 
a number of unusual features. 

As seen from Fig. [31 those electrons from the lower 
subband (— ) that have \py\ < mX, pass three turn- 
ing points in the course of their motion in the x di- 
rection, corresponding to three solutions of the equa- 
tion dH-/dpx = 0. Two of these turning points (those 
with p = mX) correspond to simultaneous vanishing of 
both velocity components, v{x) = 0, the behavior generi- 
cally impossible in a 2DEG with the parabolic spectrum, 
H — p^ /2m. The fundamental difference in the classi- 
cal dynamics of spin-orbit-split subbands lies with the 
fact that all electron states in the lower subband with 
the same energy but different momenta stop at these 
two turning points at the same point xq, U{xq) — Ep, 
provided that \py\ < mX. Consequences of this fact for 
the anomalous behavior of the ballistic conductance have 
been discussed in Ref. [2^. 

Both singularities in the spin density vp{x) = 
and vf{x) — A, originate from the accumulation of clas- 
sical turning points from many trajectories. Let us now 
demonstrate how the stronger of the two singularities, 
vf{x) = 0, is regularized when the spin density is cal- 
culated from the solutions of the Schrodinger equation 
with the Hamiltonian ([1]). We first perform decomposi- 
tion of the wavefunction ip(x,y) into a product of the fast 
exponent e^^'^ and a slow spinor function. Components 
of the latter satisfy the Schrodinger equation for a par- 
ticle in the homogeneous field F with the effective mass 
TTieff = (mX/pxYm. Corresponding solutions are the well 
known Airy functions [28]. At E = Ep = we find in 
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FIG. 3: Bottom: Family of classical trajectories at i? = Ep 
for different values of py, and m = X = F. Trajectories 
for both lower (red) and upper (green) spin-orbit split sub- 
bands are shown. An example of a trajectory contributing 
to the peak in spin density at x ~ 0, Eq. is shown 
for "Pj; <g mA" (see the text). Top: The electron density 
p = ■f/'ji/'i -I- ip2''p2 for given longitudinal momentum and en- 
ergy. Here ■!/'i,2 are two eigenfunctions of the Hamiltonian 
Eq. in with Py = 0.2mX,EF = 0,F = 0.003m^A^/R, p in 
arbitrary units, x in units of h/mX. Three classical turning 
points can be seen. The interference of incoming and reflected 
waves in the upper subband causes smooth oscillations to the 
right of the inner turning point I {x > mX^/2F). At the 
other turning points, II and HI, the two kinds of oscillations 
are seen. Slow oscillations are caused by the interference of 
the incoming wave and the wave reflected at the turning point. 
Fast oscillations (wavelength ~ Ti/mX) are due to the interfer- 
ence of distant (in time) segments of the same trajectory. In- 
set: Kinetic energy (arbitrary units) T±{px) = (p±mA)^/2m 
for fixed py — 0.2mA. 
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and 



where ^ = —x m{2FX'^ /plY^^ 
tan 2a = (mX - px)/{rnX +Px)- 

As seen from Fig. [3l electron trajectories always 
bounce twice at the turning point x = 0. This leads 
to the interference of the two solutions, Eq. ((T]), with the 
opposite signs oi Px, see Fig. [3l top. Ignoring these mi- 
croscopic oscillations that occur on the scale ^ Ti/mX, 
we can write the expectation value of the z-component 
of electron spin for x <C mX^/F as follows. 



3meV d 
AttTl dx 



dz Ai^(— x/z). 



(8) 
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where x = x{2Fm/Ti'^Y/^ . In the asymptotic region x ^ 
one can average over the oscillations of 
the Airy function. This allows us to recover the singular 
behavior of the smooth spin density ([5]): (sz) ^ x^^^'^. 

The most interesting is the behavior of spin density 
close to the turning point, for |a;| ^ (2Fm/?i^)^/'^. The 
integral in Eq. ([8]) features a logarithmic singularity at 
a; = 0. This enhancement of spin-density is due to the 
electrons with ^ mX, whose wavef unctions, given by 
Eq. ((7]), oscillate rapidly and add coherently only at the 
point a; = 0. With the logarithmic accuracy the height 
of the peak of spin density is 
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Striking feature of this result is that this maximal value 
is virtually independent of the strength of spin-orbit cou- 
pling or the shape of the boundary potential (up to a 
weak logarithmic factor). 




FIG. 4: Integrated spin S{x) = s^dx (left panel) and 
spin density Sz{z) (right panel) for F — —am m 
the vicinity of the turning point x = 0. Curves for a — 
0.25,0.05,0.01,0.002 are shown in yellow, green, blue and 
red respectively. The distance x is measured in units of 
(lOOay^^'^h/mX. Integrated spin <S and spin density are 
measured in units of {W0a)^^^eVm/8h-K^ and eVm/8Xn^ , re- 
spectively. Dashed black lines stand for the semiclassical re- 
sult (|4I5|) . The approximation (|S} valid close to the classical 
turning point is shown by solid black lines. Development of 
the logarithmic maximum of the spin density at a; = for 
small values of a, Eq. ((9)1, is clearly seen in the right panel. 

Numerical calculations presented in Fig. 0] illustrate 
the emergence of the logarithmic peak when the slope of 
the boundary potential F decreases, Eq. (0). 

In Summary, we have predicted that the nonequilib- 
rium spin-Hall spin accumulation near a smooth bound- 
ary of 2DEG ballistic conductor with spin-orbit interac- 
tion develops a narrow peak at the edge, with the width 
- {h'^/mFy^^ and height given by Eq. ([9]). It is followed 
by a slow non-monotonic decay, see Fig. O This smooth 
tail of spin density persists to much larger distances, 
^ mX^ / F. The amount of spin accumulated in the peak 
(found as a maximum of the function S{x) = J Szdx) 



equals 5max « OMeV{m'^/nFy/^ > 0.04ey/A, where in 
the last inequality we utilize the fact that F < rn^X^ jfi. 
We thus conclude that the spin accumulated at the edge 
described by a semiclassical boundary potential is in- 
versely proportional to the strength of spin-orbit inter- 
action and becomes progressively larger for smoother 
slopes. This prediction can be used for experimental ob- 
servation of spin-Hall effect in realistic two-dimensional 
electron systems. 
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